Dynamical properties of a trapped dipolar Fermi gas at finite temperature 
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We investigate the dynamical properties of a trapped finite-temperature normal Fermi gas with 
dipole-dipole interaction. For the free expansion dynamics, we show that the expanded gas always 
becomes stretched along the direction of the dipole moment. In addition, we present the temper- 
ature and interaction dependences of the asymptotical aspect ratio. We further study the collapse 
dynamics of the system by suddenly increasing the dipolar interaction strength. We show that, in 
contrast to the anisotropic collapse of a dipolar Bose- Einstein condensate, a dipolar Fermi gas always 
collapses isotropically when the system becomes globally unstable. We also explore the interaction 
and temperature dependences for the frequencies of the low-lying collective excitations. 

PACS numbers: 03.75.Ss, 05.30.Fk 
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I. INTRODUCTION 

The experimental progress towards making quantum 
degenerate gases of fermionic polar molecules [l|-[3| have 
stimulated many interests in theoretical studying of 
the dipolar Fermi gases. Compared to short-range s- 
wave interaction, the long-range and anisotropic na- 
ture of the dipole-dipole interaction may give rise to 
new phenomena in degenerate Fermi gases, such as the 
anisotropic Bardeen-Cooper-Schrieffer pairing [34To| and 
the stron gly correlated quantum phases in rapidly rotat- 
ing traps [111 ll2T] and in optical lattices p~3l - Fl5j . 

For dipolar Fermi gases in normal phase, many theo- 
retical studies focus on the zero-temperature regime. In 
early stage, the momentum distribution was assumed to 
be isotropic fl6l - [l~8j . As a result, the Fock-exchange in- 
teraction was completely ignored. Miyakawa et al. [l9[ 
introduced an ellipsoidal ansatz for the phase-space dis- 
tribution function (PSDF) which allows for the deforma- 
tion in the momentum distribution. It was then found 
that Fock-exchange interaction induces a momentum- 
space deformation, which has an important impact on 
the stability of the system. The variational approach can 
be easily implemented, but it often overestimates the sta- 
bility of the system as it is unable to capture the local 
collapses. Two of us then solve the semiclassical theory 
by numerically minimizing the total energy of the sys- 
tem [20J. In addition, Chan et al. 2l| have developed an 
analytical Fermi liquid description for the dipolar Fermi 
gas. It was also predicated that this system may support 
biaxial nematic phases [22l [23| and an interaction-driven 
quantum phase transition from a paraelectric to a ferro- 
electric quantum gas [23 |. 

As to the dynamical properties of the zero-temperature 
dipolar Fermi gases, Sogo et al. [25| have investigated the 
collective excitations and free expansions based on the el- 
lipsoidal ansatz for th e eq uilibrium distribution function. 
Lima and Pelster [26|, [27| have presented a detailed the- 
ory on the variational time-dependent Hartree-Fock ap- 
proach, which was subsequently used to study the static 



and dynamical properties of a dipolar Fermi gas in the 
hydrodynamic regime. Ronen and Bohn [28j have also in- 
vestigated the zero sound propagation in a homogeneous 
dipolar Fermi gas. 

For hnite temperature systems, two of us have studied 
the properties of an equilibrium dipolar gas through a 
full numerical calculation [29[. Baillie and Blakie have 
investigated the first- and second-order correlation prop- 
erties [30l | . In high-temperature (non-degenerate) regime, 
Endo et al. [3lj | have also introduced a variational ansatz 
which can describe the deformations in both real- and 
moment-space. So far, only the equilibrium state proper- 
ties have been explored in the finite-temperature regime. 

Experimentally, the detections of the dipolar effects 
in Bose-Einstein condensates often rely on their dy- 



namical behaviors, for instance, the free expansion 32 1 
and the low-lying collective excitations [33[. For zero- 
temperature dipolar Fermi gases, the dynamical proper- 
ties were theoretically investigated in Ref. [H, Ull based 
on the variational calculations. In this paper, we study 
the dynamical properties of a hnitc-temperature dipo- 
lar Fermi gas, to extend the works in Refs. [25|, [29| . 
Using semiclassical theory, our system is described by 
the PSDF, whose dynamical behavior is governed by the 
Boltzmann-Vlasov equation. In order to make the nu- 
merical calculation manageable while still maintain the 
necessary accuracy, we will employ a hybrid approach. 
Namely, the equilibrium distribution function is obtained 
via a full numerical calculation; while the dynamical 
equation is solved by making use of the scaling ansatz. 
Our particular interests are focused on three types of the 
dynamical behaviors: free expansion, collapse dynamics, 
and low-lying collective excitations. We will show how 
these dynamical properties depend on temperature and 
dipolar interaction. 

The content of this paper is organized as follows. In 
Sec. [Til we present our model Hamiltonian for a trapped 
single-component dipolar Fermi gas and outline the semi- 
classical description for the equilibrium state of the sys- 
tem at finite temperature. In Sec. IIII1 we introduce the 
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Boltzmann-Vlasov equation and its scaling solution. We 
then derive a set of equations for the scaling parame- 
ters which describe dynamical behaviors of the system. 
The numerical results on free expansion, collapse dynam- 
ics, and low-lying collective excitations are presented in 
Sec. IIVI Finally, we conclude in Sec. [V] 



II. PHASE-SPACE DISTRIBUTION FUNCTION 
OF AN EQUILIBRIUM STATE 

We consider a system of N single-component fermionic 
polar molecules trapped in a harmonic potential 

U ho (r) = -M {uo 2 x x 2 + u 2 y y 2 + u> 2 z z 2 ) , 

where M is the mass of the molecule and uj n (j) — x, y, z) 
are the trap frequencies. For simplicity, the trapping po- 
tential is assumed to be axially symmetric with lo x — 
Lu y = oj±. The shape of the trap is then characterized 
by the aspect ratio A = lo z /lo^_. We assume that each 
molecule possesses a permanent electric dipole moment 
d which is polarized along the positive z-axis by an ex- 
ternal electric field, such that the inter-particle dipolar 
interaction potential becomes 



V d \ 



Cd 



y 



2z 2 



(x 2 +y 2 + z 2 ) 5 / 2 ' 
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where the dipolar interaction strength is characterized 
by Cd = d 2 j (47T£o) with e being the permittivity of free 
space. Since the s-wave scattering length vanishes for 
spin polarized fermions, the Hamiltonian for the system 
under consideration takes the form 



H 



N r 

£ 

i=l 



2M 



1 N 

- ^(r i -r J -).(2) 



Within the framework of semiclassical theory, the sys- 
tem is described by the PSDF / (r, k) , which, for an equi- 
librium state, satisfies the Fermi-Dirac statistics 



/(r,k) 



exp 



e (r, k) - (j. 
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(3) 



where fc^ is the Boltzmann constant, T is the tempera- 
ture, and /i is the chemical potential introduced to con- 
serve the total number of particles 



N = (2tt)~ 



drdk/(r,k), 



and 



r(r,k) 



2u2 



ft 2 k 

~2M 



+ U e g(v, k) 



(4) 



(5) 



is the quasi-particle dispersion relation. The effective 
potential U c g contains the contributions from the exter- 
nal trapping potential and the mean field induced by the 



inter-particle dipolar interaction 



t/ cff (r, k) = U ho (r) + J drV d (r - r') n(r') 



dk' 



^(k-k')/(r,k'), (6) 



where n (r) = (2ir) 3 J dkf (r, k) is the real-space den- 
sity and 

f> n v 4tt k 2 x + k 2 y - 2k 2 

V d (k) = ~c d - kl + kl + ^ 

is the Fourier transform of Vd(r). The second and third 
terms on the right-hand-side of Eq. ([6]) originate, respec- 
tively, from the Hartree-direct and Fock-exchange inter- 
actions. 

Equations ([3|)-((6]) form a closed system of equations 
which can be solved numerically through an iterative pro- 
cedure to obtain an equilibrium PSDF. 



III. BOLTZMANN-VLASOV EQUATION AND 
SCALING ANSATZ 



The dynamical behavior of the system is described by 
the time-dependent PSDF /(r, k, t), which, in the col- 
lisionless regime, satisfies the Boltzmann-Vlasov kinetic 
equation [25|, |34[ 

1 d d 
- ^U eS (r^,t)-—f(r,k,t)=0, (7) 

Here, the effective potential J7 e ff(r, k, t) is also deter- 
mined by Eq. © with n(r) and /(r, k) being replaced 
by the time-dependent ones. 

Directly evolving the Boltzmann-Vlasov equation in 
phase space requires tremendous numerical efforts. To 
simplify the calculation, we make use of the scaling 
ansatz which assumes that Eq. ([7]) admits a solution 
of the form 



/(r,k,t) = /o(R(t),K(t)), 



(8) 



where /o(r, k) = /(r, k, t = 0) represents the equilibrium 
PSDF and 



R v (t) 



b v (ty 



Mt 



K v (t) = b v (t)kr, - —b n {t)r lv 



(9) 
(10) 



with b v being the time-dependent scaling parameters. 
Previously, this scaling ansatz has been widely adopted to 
study the dynamical properties of the quantum gase s [18. 
miSHii]- Alternatively, Lima and Pelster [H,]^ have 
presented a variational time-dependent Hartree-Fock the- 
ory for dipolar Fermi gas by employing a common- 
phase approximation. We show that the scaling ansatz, 
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Eqs. (l8|)- (fTU|) . can be derived by adopting a harmonic 
ansatz for the common phase in App. 

After introducing the scaling ansatz, the time depen- 
dence of the system is completely characterized by three 
scaling parameters b n (t). Following the standard proce- 
dure, it can be shown that the scaling parameters obey 
the following coupled dynamical equations 



h2 ( K v) , 77(b) + T/ (b) 



0, (11) 



where (K*) Q = (27r)~ 3 J dRdKK 2 f (R, K) and 

(R 2 ) = (2?r)~ 3 J" dKdKR 2 f (R, K) represent, respec- 
tively, the average sizes of the equilibrium cloud in mo- 
mentum and real spaces. The second and third terms on 
the left-hand-side of Eq. (fTTj) originate from the external 
trapping potential and kinetic terms, respectively. The 
contribution from Hartree-direct interaction is 

T " ^ = 7^3 W v (b; P) no (P) h (-P) , (12) 

z J (27T) 

where n (R) = (2n)~ 3 f rfK/o(R, K) and 
MP) 



h Q (P) = J"[n (R)], 

dW (b; R) 



Ri 



dR v 



with T [•] denoting the Fourier transform and 



W(b;K) 



blX 2 
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tip* 



2b 2 Z 2 



(b 2 x X 2 + b\Y 2 + b 2 z Z< 



,5/2 



(13) 



being the dipolar interaction potential under the scaling 
transformation. Finally, the term corresponding to Fock- 
exchange interaction is 



r/(b) = - 



dRdKdK' 



(27T) 

x/ (R,K')/o(R,K) 



W v (b; K - K') 



Once the equilibrium distribution /q (R, K) is known, the 
values of T,^ih) and T r f (h) can be calculated through 
numerical integrations for any given b. 

One can also study the collective excitations of a 
trapped dipolar Fermi gas by linearizing the dynami- 
cal equations (| 1 1 [) around an equilibrium distribution 
(b = 1). To this end, we assume that 



(15) 



with b^ being small quantities. Substituting Eqs. (|15[) 
into Eqs. (jlip and keeping the linear terms in bj, we find 
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FIG. 1: (Color online) Schematical plot of the three shape 
oscillation modes. 



where 



S, 



drdr'vj W) (r - r') n (r') n (r) 



drdkdk' 
1^? 



-r 



(W) 



(k-k')/o(r,k')/o(r,k) 



with V r 



(W) 



V, 



(vv') 



(r) = rrf-^cPVd (r) /(dr^dr^) and 



d ( k )=-^ 



(»?»?') 



(r) 



We then look for the stationary solutions of the form 

b v (t) = 6„ (0) exp (-iffi) . (17) 

Submitting Eq. (IT7|) into the differential equations (|16l) . 
we obtain a system of linear equations whose non-trivial 
solutions are determined by the characteristic equation 



-r> 2 



U,y X 

®>zx 



n 2 



'■■xy 



l yy 



<-zy 



&xz 

-tt 2 + a z 



0, (18) 



where the matrix elements are defined as 



3T„ 



2M(R 2 ) Q 



(19) 



with % = %\ H) {1) 



Tr) , (1). Utilizing the axial sym- 



(14) metry of the system, one finds a x 



'"yx 5 



and a zx = a zy , which allows us to reduce the 
number of independent matrix elements. 

Solving the characteristic equation analytically, we find 
three eigenfrequencies 



^1,2 — 2 \ ax 
o 2 — 



*xy 



a zz ± VA 



yx 



where A 



a zz f + 8a xz a z 



(20) 
(21) 

The corre- 



sponding unnormalized eigenvectors are 



ui, 2 = («i,2, «i,2, 1) T and u 3 = (-1, 1, 0) T 

with V\ t 2 = (a>xx + &xy — a zz ± y/A) /(4o ot ). In Fig. [B 
we graphically illustrate three shape oscillation modes. 
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In a spherically symmetric system, a collective oscilla- 
tion mode can be uniquely characterized by the princi- 
pal (n), azimuthal (I), and magnetic (m) quantum num- 
bers. Since our system only possesses an axial sym- 
metry, n and I are no longer good quantum numbers. 
Therefore, modes 1 and 2 result from the coupling of 
the monopole (n = 1,1 = 0,rn = 0) and quadrupole 
(n = 0, 1 — 2, m = 0) modes; While mode 3 represents 
the m — 2 quadrupole mode. 

In App.|Bl we present an alternative calculation for the 
frequencies of the three shape oscillation modes based on 
the sum-rule approach. We analytically show that these 
two approaches generate the same oscillation frequencies. 



IV. RESULTS 

To present our results, we introduce a set of dimen- 
sionless units based on the geometric average of the trap 
frequencies uj = (w^) 1 / 3 and the harmonic oscillator 
length a = y/h/(MQ): 7V 1/6 a for length, N^a' 1 for 
wave vector, and iV 1 / 3 ^; for energy. Under these choices, 
the dipolar interaction strength is characterized by a di- 
mensionless quantity 



D 



JVV6 



hwa 3. 



(22) 



In a recent experiment Q , a gas ofiV~3.9xl0 4 ground 
state KRb molecules is realized. The transverse and ax- 
ial frequencies of trapping potential are uj±_ = (2tt) 32 
Hz and u> z = (2tt) 195 Hz, respectively, resulting in a 
pancake-shaped trap potential with aspect ratio A ~ 6.1. 
From these experimental parameters, one finds the di- 
mensionless dipolar interaction strength to be 



D ~ 9.62d 2 



(23) 



where d is the electric dipole moment in units of Debye. 
For rotating polar molecules, the value of d is tunable 
by an external electric field, however it is limited by the 
permanent dipole moment (0.57 Debye for the ground 
state KRb molecule). Therefore, the maximal value of 
D can be realized in this system is around 3. Finally, 
the temperature of the gas realized in the experiment is 
T ~ 220 nK ~ 1.3T F , where T F = (6Ny/ 3 tud fk B is the 
Fermi temperature. 

Now, our system is completely specified by three pa- 
rameters: the trap aspect ratio A, the dimensionless in- 
teraction parameter D, and the temperature T (29| . In 
the following, we investigate the dynamical properties of 
a dipolar Fermi gas, including free expansion, collapse 
dynamics, and low-lying collection excitations. To this 
end, we first obtain an equilibrium PSDF fo (R, K) via 
the iterative procedure described in Sec. HU We then 
numerically evolve Eqs. (|lip with the initial conditions 
6^(0) = 1 and b v (0) — 0. Based on the scaling ansatz, 
the PSDF /(r, k, t) can then be found from Eqs. (|g|)- (flH . 
which allows us to calculate various physical quantities. 
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FIG. 2: (Color online) Time dependences of the rms cloud 
widths (a) and the aspect ratios of the distributions in real 
and momentum spaces (b) for iV = 3.9 x 10 4 , A = 6.1, ui = 
(2tt) 58.4 Hz, D = 2, and T = 0.55T F . 



A. Free expansion 

Let us first consider the free expansion of an initially 
trapped dipolar Fermi gas. As a diagnostic tool, the time- 
of-flight imaging has been used extensively in cold atom 
physics. The expanded cloud directly reflects the momen- 
tum distribution of the system, which bears the signature 
of the underlying dipolar interaction. In fact, the dipolar 
effect in Cr condensates was first detected by comparing 
the aspect ratios of the expand clouds [32| • Particularly, 
for a gas of KRb molecules, D. Wang et al. [H have exper- 
imentally realized a scheme for direct absorption imaging 
of an ultracold polar molecular gas at arbitrary external 
electric or magnetic field. 

The expansion dynamics can be studied by removing 
the restoring force terms u>^b v in Eqs. (|11|) , which corre- 
sponds to turning off the external trapping potential. In 
terms of the scaling parameters b v , the root-mean-square 
(rms) cloud widths are 



^=b n (t)yf{^~ . 



(24) 



In particular, if one switches off the inter-particle dipolar 
interaction and lets the cloud expand ballistically, equa- 
tions (fTTj) can be solved analytically to yield 



6<°> (*) = ^l + M-W^V^. (25) 
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limit. This can be easily understood by noting that the 
shape of the expanding cloud is essentially determined by 
the momentum distribution when the dipolar interaction 
is negligible. In particular, for ballistic expansion during 
which the inter-particle interaction is absent, it can be 
easily verified from Eq. (1251) that = K m (0). 

In Fig. [3J we plot the dipolar interaction strength and 
temperature dependences of the asymptotic aspect ratio 
Koo for various initial trap geometries. As a compari- 
son, the initial momentum-space deformation K m (0) is 
also plotted. A general observation is that Koo < 1 un- 
der all situations, indicating that the expanding cloud 
eventually becomes cigar-shaped. In addition, increas- 
ing interaction strength or lowering temperature results 
in larger anisotropy of the expanding cloud. During free 
expansion, both Hartree-direct and Fock-exchange inter- 
action energies are eventually converted into the kinetic 
energy. As the direct dipolar interaction tends to stretch 
the expanding cloud along z-axis, one finds that Koo is 
significantly smaller than Kk(0) in a highly anisotropic 
trap. However, in a spherical trap, where the direct dipo- 
lar interaction is small, the discrepancy between and 
K m (0) is negligible. 



FIG. 3: (Color online) Dipolar interaction strength (a) and 
temperature (b) dependences of the asymptotic aspect ratio 
Koo of the cloud for various trap geometries. The solid and 
dashed lines correspond to, respectively, the free and ballistic 
expansions. Other parameters are T/Tf = 0.55 in (a) and 
D = 1 in (b). 

We point out that, for rotating polar molecules, one can 
easily switches off the dipolar interaction by removing the 
external electric field. In Fig.[2ja), we present the typical 
results for the rms cloud widths as functions of time for a 
dipolar Fermi gas initially trapped in a pancake-shaped 
potential using experimental parameters. It can be seen 
that, similar to the zero-temperature case [l8j], the real- 
space density distribution eventually becomes stretched 
along z-axis during the expansion due to the anisotropic 
dipolar interaction. As we shall show, this conclusion is 
independent of the trap geometry. 

The axial symmetry of the system implies b x — b y = 
which allows us to characterize the time-dependent 
deformations of the distributions using 

Kr(t) = V(x 2 )/(z 2 ) and n m (t) = VWJW) 

in the real and momentum spaces, respectively. For an 
initial equilibrium distribution, K r (0) strongly depends 
on the trap geometry (In fact, for a trapped ideal Fermi 
gas, it can be easily verified that, independent of the 
temperature, the initial real space aspect ratio is exactly 
the trap aspect ratio A); while, as a result of the Fock- 
exchange interaction, the initial momentum distribution 
is always stretched along z-axis [« m (0) < 1]. For the dy- 
namical behaviors, as shown in Fig. [U[b), both n r (t) and 
K m (t) approach the same asymptotic value Kqo at large t 



B. Collapse dynamics 

Now, we turn to study the collapse dynamics of a 
trapped dipolar Fermi gas. A Fermi gas becomes unsta- 
ble when the degenerate pressure is unable to balance the 
attractive inter-particle interaction. The collapses of an 
ultracold Fermi gas was demonstrated in boson-fermion 
mixtures [40j , in which the instability of the Fermi gas is 
induced by the attractive interaction between the bosons 
and fermions. However, for a two-component Fermi gas 
with interspecies s-wave interaction, collapse has not yet 
been achieved experimentally. In fact, as the total energy 
remains positive [4l|, [42[ , the system is stable even in the 
unitary regime where the scattering length is negative in- 
finity. Only at very high density, when the inter-particle 
spacing becomes comparable to the effective range of the 
interaction, the system can in principle becomes unsta- 
ble m. 

On the other hand, with the long-range and partially 
attractive inter-particle interaction, the dipolar Fermi 
gases may provide a prospective platform to study the 
collapse dynamics of the fermionic gases. Calculations 
within mean-field have shown that, for given A and T, 
there exists a critical dipolar interaction strength D* such 
that the system becomes unstable when D > D* 29] . 

To reveal the dynamical process of a collapse, we first 
prepare an initially stable state under a given set of con- 
trol parameters: A, T, and D = Dj. At time t = + , 
the dipolar interaction strength D is suddenly increased 
to Df (> D*). We then numerically evolve Eqs. (fTT|) to 
simulate the dynamics of the system. In terms of b^, col- 
lapses occur when at least one of the scaling parameters 
approaches zero. Additionally, we say that a collapse is 
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FIG. 4: (Color online) (a) Typical dynamical behaviors of b± 
(solid lines) and b z (dashed lines) for A = 1, T = 0.287V, 
Di = 2, and various -D/'s. (b) E(b,0) /Ei,K as a function of b 
for 7,; = —0.3 and various Df /ZVs. The horizontal dotted line 
denotes the total energy E(l,0)/Ei, K for D//A = 1.3252. 



isotropic if all scaling parameters go to zero simultane- 
ously; otherwise, it is anisotropic. 

Figure 0Ja) plots the typical dynamical behaviors of 
the scaling parameters for A = 1, T = 0.28Tp, and af- 
ter the dipolar interaction strength is suddenly increased 
from Di — 2 to various -D/'s. As can be seen, the sys- 
tem becomes unstable and starts to collapse isotropically 



only when Df is larger than a threshold value D 



f 



10. 



Surprisingly, we find that DJ is significantly larger than 
D* = 2.04 for the equilibrium state. To find an expla- 
nation, we assume that b±(t) ~ b z (t) = b(t) as shown 
in Fig. |U[a). As a result, the total energy for t > + 
becomes 



E(b,b) 



E 



b 2 



' ^ - ' DfElJ >r^M Hi 



E hP b 2 



Di 6 3 



where E i>K = (2tt)- 3 / drdk (h 2 k 2 /2M) / (r,k), E uP 
(2^)~ 3 / drdkUh (r)/ (r, k), and 

E u = i J drdr'n (r)n {r')V d (r-r') 



drdkdk' 

(2tt) 6 



/o(r,k)/ (r,k')Ki(k-k') 



are, respectively, the kinetic, potential, and interaction 
energies of the equilibrium state (with dipolar interaction 
strength Di). 

Let us first consider the initial state by taking Df = Di 
in Eq. (f26|) . Consequently, the total energy is denoted as 
Ei. The stationary condition, ^j^\ b=1 j,_ = 0, yields 



2E, 



K 



2E< 



o, 



which is exactly the Virial theorem. Furthermore, the 

f) 2 F 

stability condition for the initial state, ' 
requires 



db 2 1 6 



=1,6=0 > °- 



7i 



Eij » _ 8 



indicating that the ratio of the interaction energy to the 
kinetic energy must be larger than a critical value 7* for 
the initial state. 

After the dipolar interaction strength is switched to 
D f , the stability of the system can then be analyzed by 
examining the 6-dependence of E(b, b) at t = + , i.e., 



£7(6,0) 



E, 



i,K 



= 7^+1 



< 27 > 



Apparently, the system is always stable if 7^ > 0. There- 
fore, in order to induce a collapse, one must have 7; < 
for the initial state. Moreover, to find the threshold D^, 
we plot, in Fig. IDJb), the typical behaviors of E(b, 0) cor- 
responding to 7^ = —0.3 and various D/'s. As can be 
seen, when Df/Di is not very large, £7(6,0) has a local 
maximum at 6 = 6 



the conditions 



-'max; 
dEjbfi) 
db 



with 6 max being determined by 



and 



b=b„ 



d 2 E(b,0) 
db 2 



< 



b=b„ 



0. However, this local maximum vanishes for sufficiently 
large Df/Di, under which the system becomes unstable. 
From above analysis, it becomes clear that the equation 



£7(1,0) = E{b max ,0) 



(28) 



should be satisfied when the threshold Dj is reached. 

-0.0692. 



For the example in Fig. QJa), we have 7i 
A threshold D'j = 9.68 can then be determined from 
Eq. (|28p. in agreement with our numerical finding. 

To quantitatively show that D*f is larger than D*, we 
need to find the relation between 7$ and Di numerically. 
Here, for simplicity, we present a qualitative argument. 
To this end, we consider two limit cases with 7^ — > 0~ and 
7*, for which the threshold dipolar interaction strengths 
(26) are, respectively, Dj — 00 and D* . Since D* f > D* in 
both cases, as a natural generalization, it should be held 
for any 7 e (7*, 0). 

Figure [SJa) shows the typical dynamical behaviors of 
b n for A = 0.1, T = 0.552>, A = 2, and D/'s. As 
can be seen, there also exists a threshold Dj- in a cigar- 
shaped trap. In addition, the system collapses faster with 
a larger Df. In Fig. EJb), for the fixed Di = 2 and 
Df = 10, we compare the collapse dynamics by varying 
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FIG. 5: (Color online) Typical dynamical behaviors of b± 
(solid lines) and b z (dashed lines) in a cigar-shaped trap 
(A = 0.1) after the dipolar interaction strength is suddenly 
increased from Di = 2 to Df. (a) Df dependence for a given 
temperature T — 0.55Tf. (b) Temperature dependence for 
given Df — 10. 
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FIG. 6: (Color online) Dipolar interaction strength (a) and 
temperature (b) dependences of the shape oscillation frequen- 
cies for a cigar-shaped trap with A = 0.1. Modes 1, 2, and 3 
are labeled with solid, dashed, and dash-dotted lines, respec- 
tively. 



the temperature of the system. For the high temperature 
case (T = 0.83T F ), the equilibrium PSDF of the trapped 
dipolar Fermi gas becomes more resembling of that of a 
trapped ideal Fermi gas, such that the dipolar interaction 
energy is negligible small under given Di. Consequently, 
the system remains dynamically stable. However, col- 
lapses are realized when the temperature is lowered. 

An important feature revealed in Fig. [S] is that, when- 
ever the system becomes dynamically unstable in a cigar- 
shaped trap, the scaling parameters b± and b z always go 
to zero simultaneously, suggesting that a dipolar Fermi 
gas always undergoes isotropic collapses. This feature 
is in striking difference with the anisotropic collapse of 
a dipolar condensate [II], and it is caused by the Fock- 
exchange interaction. In fact, if we artificially remove 
Tf (b) term when evolving Eqs. (|TT|) . the collapse be- 
comes anisotropic again. 

In highly pancake-shaped traps, the system remains 
dynamically stable for all control parameters we have 
tested. This can be easily understood based on our argu- 
ment for the spherical trap case. Indeed, we find that the 
total interaction energy for an equilibrium state is always 
positive. Even though the variational calculation shows 
that the interaction energy may eventually becomes neg- 
ative for sufficiently large Z^, the full numerical calcula- 
tion, on the other hand, indicates that the system col- 
lapses locally under such dipolar interaction strength. 



C. Collective excitations 

Finally, we consider the low-lying collective excita- 
tions of a trapped dipolar Fermi gas, which also re- 
quires knowing the equilibrium distribution function for 
a given set of control parameters. In the non-interacting 
limit, these oscillation frequencies reduce to 0,^1 = 

\/2uj_^l + A 2 ± | A 2 — 1| and tt^ = 2u±, which are in- 
dependent of the temperature. After the interaction is 
switched on, the collective excitation frequencies will be 
shifted with respect to ■ For the same D value, the 
interaction induced frequency shifts in highly anisotropic 
traps are usually larger than those in a spherical trap, 
due to the strong direct dipolar interaction. Therefore, 
we will concentrate on the low-lying collective excitations 
in cigar- and pancake-shaped traps. 

Figure |6]Ja) shows the dipolar interaction dependence 
of the shape oscillation frequencies for A = 0.1 and 
T = 0.055Tp. For relatively weak dipolar interaction, 
the frequencies of all three modes are slightly shifted up- 
wards, whereas, the frequency of the mode 1 (monopole) 
starts to go down for D > 1.2. Close to the stability 
boundary, the value of J7i drops significantly. The fact 
that the monopole mode goes soft in strong dipolar inter- 
action regime is in agreement with the variational calcu- 
lation at zero-temperature limit [2o| . Moreover, it is also 
consistent with the isotropic collapse in a cigar-shaped 
trap discussed previously. We remark that the softening 
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of the monopole mode sensitively depends on tempera- 
ture of the system. For instance, when the temperature 
is increased to T = 0.55Tp, the frequencies of the three 
modes shift all the way upwards. This can be under- 
stood by examining the equilibrium distribution function 
at the vicinity of D*, from which one may identify the 
type of the instability when the system becomes unsta- 
ble. In fact, for the low temperature case (T = 0.0552f), 
the volume of the equilibrium gas goes to zero when D 
approaches D* , indicating that the system undergoes a 
global collapse. However, for the high temperature case 
(T = 0.55Tp), local collapse is induced at the onset of 
the instability, which is not described by the three shape 
oscillation modes discussed here. 

In Fig. |6jb), we plot fli as functions of T for A = 0.1 
and D = 1. Among the three shape oscillation modes, 
the frequency of the mode 3 has the largest deviation 
from that of an ideal gas, for which the dipolar inter- 
action induced frequency shift can be as high as 10% 
at the low temperature limit. When the temperature is 
increased, the frequencies of all shape oscillation modes 
shift downward and asymptotically approach those of a 
non-interacting gas. 

For pancake-shaped traps (A = 6.1), Figure [7ja) 
shows the D dependences of the shape oscillation fre- 
quencies. Taking into account the rescaling frequencies 
£lf\ the highest mode has a breathing geometry (mode 
1), whereas the lowest one is the m — 2 quadrupole 
mode. Even though the frequencies of all shape oscilla- 
tion modes shift downwards, we do not find the frequency 
of any particular mode drops significantly close to the in- 
stability, which is the manifestation of the local collapses 



for the given trap geometry. The temperature depen- 
dences of the mode frequencies are presented in Fig.[7Jb). 
As one increases the temperature, f2j increase monoton- 
ically to approach the frequencies corresponding to an 
ideal gas, showing an opposite tendency compared to the 
cigar-shaped trap. Again, the m — 2 quadrupole mode 
has the largest frequency shift with respect to an ideal 
gas. 



V. CONCLUSIONS 

In this paper, we have explored the dynamical proper- 
ties of a trapped dipolar Fermi gas at finite temperature. 
For free expansion, we find the expanded cloud always 
becomes stretched along the direction of the dipole mo- 
ment. We further explore the temperature and inter- 
action strength dependences of the asymptotical aspect 
ratio of the expanded cloud. We have also studied the 
collapse dynamics by suddenly increasing the dipolar in- 
teraction strength. In contrast to the anisotropic col- 
lapse of a dipolar condensate, we find that dipolar Fermi 
gases always collapse isotropically. Finally, we have in- 
vestigated the low-energy shape oscillations of a trapped 
dipolar Fermi gas. It is shown that, in a cigar-shaped 
trap, the monopole mode goes soft close to the instabil- 
ity, which is consistent with the isotropic collapse of the 
dipolar Fermi gas. In addition, among three shape oscil- 
lation modes, the m = 2 quadrupole mode always has the 
largest interaction induced frequency shift with respect 
to that of an ideal gas. 
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Appendix A: Scaling ansatz and common phase 
approximation 

It is well-known that the phase-space distribution func- 
tion /(r, k) is related to the single-particle reduced den- 
sity matrix as 



(Al) 



where pi (r, r') = Tr p^ (r') ip (r) 



with p being the 



density operator and tp the annihilation operator for 
fcrmionic field. To obtain an explicit expression for pi, 
we assume that {<fi a } a set of single-particle orbitals, ob- 
tained by self-consistently solving the Hartree-Fock equa- 
tions. For an equilibrium state, the single-particle re- 
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duced density matrix then becomes 

<t>a (r)P a (r') 



Pi (r,r') 



E 



(A2) 



where e a is the energy of the a-th single-particle orbital. 

We remark that, from the above equation, one may 
derive the Fermi-Dirac distribution for a trapped Fermi 
gas, Eq. by adopting the local density approxima- 
tion. To this end, we assume that, in the vicinity of a 
given space point, the system is treated as a homoge- 
neous system with a spatially dependent chemical poten- 
tial (i (r) = p — Uho (r) . Thus the single-particle orbitals 
are taken to be plane waves with the subscription a be- 
ing replaced by the wave vector k. Consequently, the 
single-particle reduced density matrix becomes 



Pi (r,r') 



elk 



zk-(r— r ) 



(2nf exp [(e(^, k) - M ) fk B T] + 1 ' 



(A3) 



where the quasi-particle dispersion relation e (r, k) takes 
the same form as Eq. l[5]). The Fourier transform of 
Eq. (|A3[) is exactly the Fermi-Dirac distribution Eq. ([3]). 

To study the collective motion of a zero-temperature 
dipolar Fermi gas, Lima and Pelster employed a com- 
mon phase approximation for the single-particle orbitals 
in Ref . (2(| [27| , which assumes that 



_ JM X (r,t)/h 



4>a(r,t), 



(A4) 



where the common phase x ( r > t) represents the collec- 
tive motion of the system. The common phase is further 
assumed to adopt a harmonic ansatz J2(| 



(A5) 



Clearly, x(r, t) gives rise to a common velocity field of 
the form 



v v {r,t) = /3 v (t)r v . 



(A6) 



Without loss of generality, we may set (3 V — b n /b n . 
The velocity field, Eq. (|A6|) . can then be interpreted as 
a global dilatation of each single-particle orbital induced 
by the scaling transform Eq. @. We can now rewrite 
the wave function for the orbitals as 

<p a (r, t) = (& s & y &,)-VVJ*x(r,t)/V a (R, 0) , (A7) 

where % = 677^/67/ and the factor (bxbybz) -1 / 2 is 

introduced to maintain the normalization. Consequently, 
the time-dependent single-particle reduced density ma- 
trix becomes 

iM[x(T,t)- X (r' ,t)]/H 
Pl (T,T',t) = - — Pl (R,R,0). (A8) 

X OyO Z 



Using Eq. (jAll) . one finds the time-dependent phase- 
space distribution function 



/(r,k,i) 



dS exp 



i \b v k v — M h 1 b n r ri ^j S v 



= /(R,K,0), (A9) 

where K is defined by Eq. (flUl) . Apparently, Eq. (IA9[) is 
exactly the scaling ansatz introduced in Sec. IIIII 



Appendix B: The sum-rule approach to collective 
excitations 

Here we use the sum-rule approach to calculate the fre- 
quencies of the three shape oscillation modes [45|,[46[. We 
assume that E n and \E n ) are, respectively, the eigenen- 
ergies and eigenstates of a many-body Hamiltonian H. 
Moreover, F is an excitation operator of the system. In 
the sum-rule approach, the average excitation energies of 
the low-lying excitations, which correspond to different 
types of shape oscillations, can be obtained through 




where 



E 



E,, 



F 



E, 



Pn {En' — E n ) 



is the fc-th order moment of the strength function for 
the transition operator F, with P n being the probability 
distribution of the eigenstate at given temperature T. 
Alternatively, we may express mi and 7713 as 



mi 



m 3 



H,F 
H 



}■ 



H. 



H, F 



where p = J2 n Pn\E n ) (E n \ is the density operator of the 
system. 

In order to obtain the frequencies of the shape oscilla- 
tions, we consider excitation operators 



F (n, I, m) = rt n+l Yim{0i, <k) , 

i=l 

where (ri,8i,(f>i) is the spherical coordinate of the z-th 
particle. Since modes 1 and 2 result from the coupling 
of the monopole (n. = 1, 1 = 0, m = 0) and quadrupole 
(n = 0, 1 = 2, m = 0) modes, we construct an excitation 
operator which is a linear combination of -F(1,0, 0) and 
F(0,2,0),i.e., 

N 

F a = [sin a (xl- + yf) + cos azf] , 
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where a G [0,7r) is a parameter to be determined. Af- 
ter a straightforward calculation, we obtain the average 
excitation frequency as a function of a, 



Sip (a) = 



i £1 cos 2 a + £2 sin a cos a + £3 sin 2 a 



cos 2 a + £ sin 2 a 



(Bl) 



where & = 4w 2 A 4 / 3 + (3% + S ZZ )/(2M (z 2 ) Q ), £ 2 = 
2& z /(M(z 2 ) ), £3 = 8w 2 A~ 2 / 3 K 2 (0) + (3% + + 
S xy )/(M(z 2 j Q ),a,ndC = 2K 2 r (0). 

To determine a, we consider two general excitation op- 
erators, F\ and F2, with average excitation frequencies 
£lp i < £lp 2 . In addition, we assume that the modes ex- 
cited by F\ and F2 are orthogonal. It can then be seen 
that, for an operator F constructed by superimposing F% 
and F2, we always have Sip <Slp < flp • Therefore, for 
eigen- excitation modes, such as those described by Eq. 
(|22[) . the values of a should either maximize or minimize 



Sip 



(a), which yields 



COS OL\ 2 



£f + - 6 ± V«2 2 + (CCi - 6) 



(B2) 



Submitting Eq. (|B2I) into (|B1|) . we obtain the frequencies 
of the shape oscillation modes 1 and 2, which take exactly 
the same form as those in Eq. (l20l) . 



The shape oscillation mode 3 can be directly excited 
by the transition operators F(0, 2, ±2). The average ex- 
citation frequencies are exactly Eq. (f2"Tj) and they are 
degenerate for m = ±2 modes due to the axial symmetry 
of the system. 
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